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We discuss the links between Generalized Parton Distributions (GPDs) and elastic nu- 
cleon form factors. These links, in the form of sum rules, represent powerful constraints on 
parametrizations of GPDs. A Regge parametrization for GPDs at small momentum trans- 
fer, is extended to the large momentum transfer region and it is found to describe the basic 
features of proton and neutron electromagnetic form factor data. This parametrization is 
used to estimate the quark contribution to the nucleon spin. 
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I. INTRODUCTION 

Generalized parton distributions (GPDs) are universal non-perturbative objects entering 

the description of hard exclusive electroproduction processes (see Refs. for reviews and 

references). These GPDs, which are defined for each quark flavor (u, d, s), parametrize nonforward 
matrix elements of lightcone operators. They depend upon the longitudinal momentum fractions 
of the initial and final quarks and upon the overall momentum transfer t to the nucleon. When the 
momentum fractions x + x — ^ of initial and final quarks are different being the longitudinal 
momentum asymmetry, or skewness), one accesses quark momentum correlations in the nucleon. 
Furthermore, if one of the quark momentum fractions is negative, GPDs reflect an antiquark 
contribution, and consequently one can investigate qq configurations in the nucleon. Therefore, 
these functions contain a wealth of new nucleon structure information, generalizing that obtained 
from inclusive deep inelastic scattering. 

In hard exclusive processes, such as deeply virtual Compton scattering, GPDs enter in most 
observables through convolution integrals. Hence, to access GPDs, the most realistic strategy to 
date seems through judicial parametrizations. Building self-consistent models of GPDs is, however, 
a rather difficult problem, because one needs to satisfy many physical principles and constraints 
which should be obeyed by GPDs. They include spectral properties, polynomiality condition, 
positivity, relations to parton densities and form factors 

mm 

In this paper, we elaborate on the i-dependence of the ^ = generalized parton distributions, 
and its interplay with the x-dependence. This subject has attracted a considerable interest. In 
particular, it has been shown j^, [lO, 11| that by a Fourier transform of the f-dependence of GPDs, 



it is conceivable to access the spatial distribution of partons in the transverse plane, and to provide 
a 3-dimensional picture of the nucleon [l^ . [l^ . The f-dependence of moments of GPDs has also 
become amenable to lattice QCD calculations recently. As the lattice calculations mature 
further, they may eventually provide additional constraints on moments of generalized parton 
distributions. Phenomenological estimates of the ^-dependence and t-dependent parametrizations 
of GPDs have already been discussed in Refs. 10, 13, 0, lis, ll9[ , and more recently, in 

Ref. [20]. Some results of the present paper were reported in Refs. 2l|, l22l |. 

We give here several parametrizations of the f-dependence of the GPDs, both at small and large 
values of —t (with t < 0, i.e. in the spacelike region). We start in Section [TI] by reviewing the 
relevant sum rules which link GPDs to form factors. Subsequently, we discuss in Section IIIII a 
Gaussian ansatz for the t-dependence of GPDs (at large —t) which has been introduced and used 



in Refs. jl5l. Ilq. Such a Gaussian ansatz, however, is not able to describe the small —t behavior 
of GPDs, and in particular gives divergent rms radii for the nucleon electromagnetic form factors. 
We therefore proceed in Section |IV] to describe a Regge parametrization 0, [21[ which provides a 
physically consistent behavior of form factors at small —t. We extend this model then in Section IVl 
to large —t so as to yield the observed power behavior of the electromagnetic form factors at large 
(spacelike) momentum transfers. We found a quite economical parametrization that allows for a 
description of both proton and neutron electromagnetic form factors with only 3 parameters: the 
universal Regge slope a' and two parameters r]u,r]d governing the a; — > 1 behavior of the splin-flip 
GPDs £^{x, t = 0),£'^{x,t = 0) relative to that of the usual parton densities u{x), d{x). We discuss 
the comparison of our results with the data in Section IVIl and use our parametrization to estimate 
the quark contribution to the nucleon spin. In Section IVIIl we discuss the positivity constraints 
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on GPDs in the impact parameter bj^ representation. To extend the region in x and where 
the positivity constraints are satisfied, we propose a model in which the parameters rj^ and rjd are 
equal. It provides (with just two parameters) almost the same quality description of the four form 
factors as the 3-parameter model. Our conclusions are presented in Section [Villi 



II. FORM FACTORS AND GPDS 

The nucleon Dirac and Pauli form factors Fi{t) and F2{t) 

F^{t) = Y.e,Ff{t) (1) 

<? 

can be calculated from the valence quark GPDs H and E through the following sum rules for their 
flavor components {q = u,d) 

r+l 



= J ^ dxH'i{x,^,t), (2) 

Flit) = j^\x E^{x,i,t) . (3) 

Since the result of the integration does not depend on the skewness ^, one can choose = in the 
previous equations. Furthermore, the integration region can be reduced to the < x < 1 interval, 
introducing the nonforward parton densities [15]: 

n'i{x,t) = H'^{x,0,t) +H''{-x,0,t), (4) 
£'^{x,t) = E'i{x,0,t) + E'i{-x,0,t), (5) 

obeying the conditions 

f\xW{x,t) = F^{t) , (6) 
Jo 

[\x£\x,t) = F!l(t) , (7) 
Jo 

that follow from the sum rules ([2]), ([3]). The 'H'^{x^t) functions also satisfy the t — > reduction 
relations 

7e{x,t = Q) = u^{x), n'^{x,t = 0) = d^{x), (8) 

connecting them with the usual valence quark densities in the proton. The t = limit of the £'^{x, i) 
distributions exists, but the "magnetic" densities £'^{x,{)) = £'^{x) cannot be directly expressed in 
terms of any known parton distribution: they contain new information about the nucleon structure. 
However, the normalization integrals 

rl 

dxS'^ix) (9) 







are constrained by the requirement that the values -Ff = 0) and F2{t = 0) are equal to the 
anomalous magnetic moments of the proton and neutron. This gives 

Ku = 2Kp + K„ PS +1.673 , (10) 
Kd = Kp + 2 K„ PS -2.033 . (11) 

For comparison, the normalization integrals for the TL^{x) = Uy{x) and Tl'^{x) = dv{x) distributions 
are given by 2 and 1 respectively, the number of u and d valence quarks in the proton. 
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III. GAUSSIAN ANSATZ 



The simplest model for the proton's Tf^^x, t) is to separate the x and t-dependencies and express 
it as the product 



W{x,t)=q,{x)Fi{t) 



(12) 



of the parton density qv{x) and the Fi{t) form factor of the proton. It trivially reproduces qv{x) 
in the forward limit and gives the correct result for Ff{t). However, such a complete factorization 
of the X and t dependencies seems rather unrealistic. In particular, the form factor formula 23 1 



oo „ n 

/ Y[d^ki^dxi'^ea'i'*p'{xi, . . . , 



3^ni fell 



xiq±, . . . , feax + (1 - Xa)qj_, ...,K 



n=l 2=1 



x^p(rci,...,x„;fci^,...,fca^,, 



(13) 



u=l 



i=l 



of the light-cone formalism is a convolution of the light cone wave functions contain- 
ing nonfactorizable combinations ki^ — Xiq±. Furthermore, the n-body Fock component 
'^p{xi, . . . ,Xn;ki^, . . . ,knj^) of the light-cone wave function usually depends on the transverse 
momenta {ki^} through the J2i kf^/xi combination involving both ki^ and the fractions x-i of the 
hadron longitudinal momentum carried by the quarks. If the dependence on this combination has 
a Gaussian form, the /cj_ integration can be performed analytically providing an example of the 
interplay between the x and t dependencies. The result of integration can be most easily illustrated 
on the simplest example of a two parton system (n = 2). In this case 



F^'^\q\) = r dxcfk±^*{x;k± + {l-x)q±)'^{x;k±) 
Jo 

Assuming the Gaussian ansatz 
^'(x; k±) ~ exp 



k', 



kl 



2xA2 2(1 - x)A2 



exp 



2a;(l -x)A2 



we obtain 



dxq^'^\x)e 



(14) 



(15) 



(16) 



where q^'^\x) has the meaning of the two-body part of the quark density q{x). This suggests the 
Gaussian (G) parametrization 15|, |2j] for the nonforward parton densities 



nUx,t)=q,{x)e^'-^^'/'^^" 



(17) 



containing a nontrivial interplay between x and t dependencies. The scale characterizes the 
average transverse momentum of the valence quarks in the nucleon. The best agreement (within 
10%) between experimental data for Ff(t) in the moderately large t region 1 GeV^ < — t < 10 GeV^ 
and calculations based on Eqs. ([I]), (l6|), (fT7|) is obtained for A^ ~ 0.7 GeV^. This value corresponds 
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to an average transverse momentum of about 300 MeV [l5i], which is close to the inverse of the 
proton size. The latter can also be estimated by calculating the mean squared radius 



"^^'P - ^ dt 



(18) 

t=o 



The Gaussian model for Ti.'^(x,t) then gives the expression 

'^i,p = 6^ dx |e„ Uy{x) + ed (it,(x)| ■ (19) 

If one assumes the standard Regge-type behavior qy{x)\x^o ~ x~^'^ of the parton densities at small 
X, the integral in (fT9]) diverges. To get a finite slope we should modify the model for 7i'^{x,t) in 
the region of small x. 

IV. SMALL T BEHAVIOR AND REGGE PARAMETRIZATION (Rl) 

The Regge picture suggests a x~°*-*^ behavior at small x or the 

7^^(x, t) = g„(x)x-("W-°(°)) (20) 

model for the nonforward densities 7i^{x,t). Assuming a linear Regge trajectory with the slope 
a , we get 

W^i(x,t) =g„(x) X-"'* . (21) 



This ansatz was already discussed in Ref. [g]. The u and d flavor components of the Dirac form 
factor are then given by 



F^{t) = f dxtx^x) e-*"''"^ , Ff{t) = f dxd^x) e-*"''-^^ . (22) 
Jo Jo 

The proton and neutron Dirac form factors follow from 

Ff(t) = euF^it) + eaFfit), (23) 
F^{t) = e^Ffit) + e^F^it). (24) 

By construction Ff (0) = 1, and Ff (0) = 0. The Dirac mean squared radii of proton and neutron 
in this model are given by 

r^p = —6 a j |eu u-^ix) + (i„(x)| Inx , (25) 

= —6a j dx |eu dv{x) + n„(x)| Inx . (26) 

Instead of the 1/x factor present in the Gaussian model, we have now a much softer logarithmic 
singularity at small x, and the integrals for rf converge. 

To calculate F2, we need an ansatz for the nonforward parton densities £'^{x, t). We assume the 
same Regge-type structure 

4i(x,t)=f'?(x)x-"'* (27) 
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as for TC'^{x,t). The next step is to model the forward magnetic densities £'^{x). The simplest idea 
is to take them proportional to the 7i'^{x) densities. Choosing 

£^{x) = yMx) and £'^{x) = Kddy{x) , (28) 
we satisfy the normalization conditions Q which, in their turn, guarantee that F2{0) = Kp, and 

As we will show in Section IVTl the Regge model Rl fits Ff(t) and i^'fl^) data for small momentum 
transfers —t ^ 0.5 GeV^. However, the suppression at larger —t in the Rl model is too strong, and 
it consequently falls considerably short of the data for —t>l GeV^. 



V. LARGE T BEHAVIOR AND MODIFIED REGGE PARAMETRIZATION (R2) 

To improve the agreement with the data at large —t, we need to modify our models. Note, that 
both the Gaussian (G) and the Regge-type model (Rl) discussed above have the structure 

7i{x,t) = qy{x) ex.p[tg{x)] , 

with g{x) ~ (1 — x)/x and g{x) ~ — Inx, respectively. Hence, at large t, the form factors are 
dominated by integration over regions where tg{x) ~ 1 or g{x) ~ 1/t ^ 0. In both cases, g{x) 
vanishes only for x — > 1, and the large-t asymptotics of Fi{t) is governed by the x ^ 1 region. 
Given g{x) ~ 1 — x as x — > 1, one derives that if Qvix) ~ (1 — x^ for x close to 1, then the form 
factors drop like l/f^^^ at large t. Experimentally, u is close to 3, thus the models G and Rl 
correspond to the ~ 1 /f^ behavior for the form factors. This seems to be in contradiction with the 
experimentally established 1/t^ behavior of Ff{t), so one may be tempted to conclude that these 
models have no chance to describe the data. A trivial but important remark is that the model 
curves for Ff{t) are more complicated functions than just a pure power behavior ~ 1/t^. In fact, 



up to lOGeV^, the Gaussian model reproduces the data for F[ within 10% 15|. For higher t, 
the Gaussian model prediction for Ff drops faster than and goes below the data. However, 
the nominal 1/t^ asymptotics is achieved only at very large values —t ~ 500 GeV^. As we show 
in Section IVTl the Regge-type model Rl result visibly underestimates the data for F^ already for 
-t ~ 1 GeV^ though one should wait tih -t ~ 100 GeV^ to see that the l/f^ behavior really settles. 
Thus, the conclusions made on the basis of asymptotic relations might be of little importance in 
the experimentally accessible region: a curve with a "wrong" large-t behaviour might be quite 
successful phenomenologically in a rather wide range of t. 

The shortcomings of the G and Rl models are more of a theoretical nature. Namely, they do not 



satisfy the Drell-Yan (DY) relation [23, l25|] between the x — > 1 behavior of the structure functions 
and the t-dependence of elastic form factors. According to DY, if the parton density behaves like 
(1 — xY, then the relevant form factor should decrease as for large t. Such a relation 

does not hold if g{x) ~ 1 — x but it holds if g{x) ^ (1 — x)^. Thus, the simplest idea is to attach 
an extra (1 — x) factor to the original g{x) functions. To preserve the Regge structure at small x 



and t we take the modified Regge ansatz R2 10|, [26 1 



n%,{x,t) = q,{x)x~'' . (29) 



The inability of the G parametrization to satisfy the DY relation may seem rather surprizing 
in view of the fact that the original derivation of the relation by Drell and Yan 23] is based on 
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the analysis of the large-g_L hmit of the general formula (|13p of which the G ansatz is a specific 
case corresponding to n = 2 and the 'if{x;k±) ~ exp [-kl/2x{l- x)\^] wave function. Note, 
that if the wave function '^{x,k±) depends on k± through the combination k'^/x + k'j_/{l — x), 
then the restriction on the x — > 1 integration region should be \k± + (1 — x)q±\'^/{l — x) ^ 
which results in the 1 — x ^ }? jq^^ constraint on the x integration. Also, from the explicit form 
of the Gaussian parametrization (jl7p . it is clear that the essential region for the Xa integration 
is 1 — X(j ~ }? l{—t) which gives the l/t'^^^ result, that differs from the canonical l/t^'^^^^/^ DY 
prediction. The resolution of this discrepancy is rather simple. In fact, in the derivation given 
by Drell and Yan, it was implied that the wave function depends on k^ through the combination 
{kP\_ + mg)/x(l — x), with being the (constituent) quark mass. Then, in the Gaussian case, after 
the A;_L-integration, one would have the structure ~ exp{ — [(1 — x)g^+m^/(l — x)]/A^]} in the x ~ 1 
region, and at large q^^ the dominant contribution comes from the region 1 — x ~ ruq/q^. This 
agrees with the argumentation of Ref. 23] , that the leading contribution to the form factor is due 
to integration over the region 1 — Xq < mq/q± where the longitudinal momentum fraction Xq of the 
active quark is close to 1 and those of the passive quarks are close to 0, so that \ka^ + (1 ~ Xa)qi\ 
and all \ki^ — Xiq_[_\ are bounded by 0(A). Integration over all ki^'^s and Xj's of passive quarks 
gives q{xa)- If q{xa) ~ (1 — XaY , then the final integration over the region Xa ~ 1 — A/(;_l gives 
F{q±) ~ 1/(Z1+^ ~ 

Turning back to the Gaussian model with zero quark mass, it is easy to realize that the factor 
(1 — x)t/xA^ in the exponent of the G parametrization may be viewed as [(1 — x)g^]^/x(l — x)A^ 
with l/x(l — x) coming from the exp[— A:^/2x(l — x)A'^] structure of the /c^-dependence of the 
wave function. As we have seen, to get the Regge-type behavior at small x, one should soften the 
1/x factor in the exponential substituting it by In x. Since the limit x^ — > 1 for the active quark 
corresponds to the Regge limit x^ ^ for the spectators, one may expect by analogy that the 
1/(1 — x) singularity is also softened after inclusion of higher Fock components. The R2 ansatz 
corresponds to substitution of the 1/(1 — x) factor by a constant. Other arguments in favor of the 



R2 model can be found in Ref. 

The correlation between the power behavior of form factors and the behavior of inclusive struc- 
ture functions W{x-q) of deeply inelastic scattering at large Bjorken variable xb is a rather popular 
subject ("inclusive-exclusive connection"). The basic idea behind the possibility of such a correla- 
tion is that, for sufficiently large x^, one approaches the exclusive single-hadron pole. The invariant 
mass W"^ = {p + q)"^ of the hadronic system produced in deep inelastic scattering is related to the 
Bjorken variable xb by 

W'^ - Mi 

l-XB = XB , (30) 

and the single-hadron contribution to cross section is given by the form factor squared multiplied by 
6{W'^ —mj^). The Bloom-Gilman duality idea ^] is that the M^^-integral of the hadron contribution 
is equal to the x-integral of the structure function Wi{x) over a duality region with fixed boundaries 
in the variable W'^. This gives a relation between the power v specifying the (1 — xY behavior 
of the structure function Wi{x) in the x — > 1 region and the power-law behavior of the squared 
elastic form factor: F'^{t) ~ (l/|t|)'^'^^. In the proton case, with usually adopted value = 3, one 
obtains a dipole behavior for the Dirac Fi{t) form factor. 

We would like to strongly emp hasize here that one should not confuse the Bloom-Gilman duality 
with the Drell- Yan relation [23|]. As we discussed above, the latter connects some integral of a 
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nonforward parton density 7i^{x^t) over the interval x > \ — \/^/^ with the first power of the 
form factor. It is worth to repeat and stress the statement: the Bloom-Gilman relation connects 
an X integral of the structure function with the square of the form factor, while the DY relation 
expresses an(other) x integral of the structure function in terms of the first power of the form factor. 
Moreover, the dominance of the region x > 1 — implied by the DY relation is a consequence 

of a specific structure of the density the interplay between its x and t dependence. As 

we have seen, the Drell-Yan relation does not work for the Gaussian model, but it holds for the 
modified Regge model R2. 

One should also realize that both relations were formulated before the QCD era, and in ab- 
solutely nonperturbative terms. Their authors did not assume that the shape of the structure 
function Fi{x) or that of the nonforward parton densities Ti.'^{x,t) are generated by perturbative 
QCD dynamics based on hard gluon exchanges. Their prescription was that knowing the x ^ 1 
behavior of the structure functions, one can use Bloom-Gilman or Drell-Yan relations to get pre- 
dictions for form factors. Both relations have a common feature: if one changes the power v in 
the (1 — x)'^ behavior of the structure function, this would result in a change of the l/(— 
power behavior of the form factor, i.e., the powers themselves are not fixed, what is fixed is the 
relation(s) between them. Accidentally, both relations give the same correlation between the two 
powers, and that is why they are confused sometimes. 

In distinction to the Bloom-Gilman and Drell-Yan relations, perturbative QCD predicts definite 
powers for the asymptotic behavior of form factors and the x — > 1 behavior of parton distributions. 
For example, it gives (as/|t|)"~^ for a spin-averaged form factor of an n-quark hadron, and it also 
predicts fixed powers a^"'~^(l — x)^"~^ for the x — > 1 behavior of its valence quark distributions 
(see [28[) 71[. The basic difference between the pQCD formulas and BG & DY relations is that a 
particular power behavior of a hadronic form factor in pQCD is not a consequence of a particular 
limiting power behavior of the respective parton distribution in the region x — > 1 . The fixed powers 
predicted by pQCD are correlated simply because of similarity of the relevant diagrams, but there 
is no causal connection between them. Also, though the powers predicted by pQCD for the nucleon 
are in agreement with BG and DY relations, it was never demonstrated that there is a fundamental 
reason behind this fact. 

Formally, the relevant powers of (1 — x), and for the proton are correlated in pQCD 
just like in the Bloom-Gilman relation. However, a direct calculation of pQCD diagrams for Wi(x) 
gives expressions which have more complicated structure than the squares of form factors (see 
e.g., [62], where the x — > 1 behavior of GPDs is also discussed). Thus, it is not clear yet if the 
Bloom-Gilman relation works in pQCD. 

With the Drell-Yan relation, the situation is simpler. The whole logic of the hard-rescattering 
pQCD mechanism is orthogonal to the Feynman-Drell-Yan approach. In the pioneering paper by 
Lepage and Brodsky |28|, it was emphasized in the Introduction of that paper that the Drell- 
Yan relation is invalid in pQCD. It was stressed, in particular, that the correlation between the 
powers of as in pQCD predictions disagrees with the Drell-Yan relation. For instance, the leading 
(1 — x)^ term in Wi{x) for the nucleon is attributed in pQCD to diagrams involving four hard gluon 
exchanges, and is accompanied hence by the factor. Integrating it over the region x > 1 — X/Q, 
one would get a contribution ~ a^/t'^ that has the same power as the pQCD prediction for 
the nucleon form factor, but has two extra powers of Og- 

Our models imply the dominating role of the Feynman-Drell-Yan mechanism for the hadronic 
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form factors, and we assume that the x ^ 1 behavior of the parton distributions is generated by 
nonperturbative dynamics. In this scenario, the observed behaviour of hadronic form factors is 
also due to the nonperturbative dynamics, and we treat as neghgible the pQCD contributions to 
the nucleon form factors, which have (as/vr)^ suppression compared to the nonperturbative terms. 

In the following estimates we take the unpolarized parton distributions at input scale iJ? = 1 
GeV2 from the MRST2002 global NNLO fit 3 as : 



= 0.262 (1 - x)3-50 [I + 3.83 x^'^ + 37.65 xj , (31) 
dy = 0.061 (1 - (l + 49.05 + 8.65 x) . (32) 

One sees that = 3.50 and = 4.03 at a scale = 1 GeV^. Hence, the asymptotic behavior 
of Ff{t) in the R2 model is 1/t^'^^, generating a slightly faster decrease than the "canonical" 1/i^. 
Again, this asymptotic limit sets in for very large t values. 

At small t, the modifications compared to the Rl model are not very significant numerically. 
The Dirac mean squared radii of proton and neutron in the R2 model are finite and given by 

rfp = — 6a J dx |eu Uy{x) + ed d^ix)^ (1 — x) Inx , (33) 

= —6a J dx |eu dy{x) + tit,(x)| {1 — x) Inx . (34) 

In case of the Pauli form factor F2, we perform the same modification of the ansatz for the 
f (x,t) densities taking 

=f'?(x)x-"'(i-^)* . (35) 

Experimentally, the proton helicity flip form factor F2{t) has a faster power fall-off at large t 
than Fi{t). Within all our models, this means that the x ~ 1 behavior of the functions £{x) and 
Ti.{x) should be different. To produce a faster decrease with t, the x ~ 1 limit of the density £'i{x) 
should have extra powers of 1 — x compared to that of W(x) (in case of the G model, such a 
modeling was originally incorporated in Ref. [l3l)- Aiming to avoid introducing too many free 
parameters, we try the simplest ansatz for £'^[x) in which we get them by just multiplying the 
valence quark distributions by an additional factor (1 — x)^', i.e., we take 

^"(3:) = ^(1 - xf-u,{x) and £\x) = ^(1 - x)^'^4(x) , (36) 

l^d 

where the normalization factors and A'^,^ 

= dx {I - xf^ Uy{x) , (37) 





Nd = I dx (1 - x)^"^ (38) 
Jo 

guarantee the conditions ([9|). The flavor components of the Pauli form factors are now given by 

F2"(t) = /^dx^(l-x)''"n„(x)x""'(i-^')* , (39) 
/■I / 



F^{t) = / dx-f (l-x)''''4(x)x-" . (40) 

Jo Nd 
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The powers rju and r]^ are to be determined from a fit to the nucleon form factor data. Note 
that the value rjg = 2 corresponds to a 1/t asymptotic behavior of the ratio F2{t) / Fl{t) at large 
t. We also tried an even simpler 2-parameter version of the R2 model, with r/„,% restricted to be 
equal to each other = 



VI. RESULTS 

In this section, we show the results for the proton and neutron electric and magnetic form factors 
based on the Regge and modified Regge parametrizations discussed in this work. In recent years, 
a lot of high accuracy data have become available for the nucleon electromagnetic form factors in 
the spacelike region, which put stringent constraints on our parametrizations of GPDs. 

The parametrization Rl of Eqs. (j21l28p depends on only one parameter : a', which can only 
be varied within a narrow range if it is to be interpreted as a slope of the Regge trajectory. The 
modified Regge parametrization R2 of Eqs. (j29|35l) depends on three parameters. Besides a', it also 
depends on rju and which govern the x — > 1 behavior of the GPD that in turn is determined 
from the behavior of -Ff /-^f at large —t. In determining these parameters, we perform a best fit 
to the Sachs electric and magnetic form factors, as they are the usual form factors extracted from 
experiment. The Sachs electric and magnetic form factors are determined from Fi and F2 as 

GE{t) = Fi{t) - TF2{t) , (41) 
GM{t) = Flit) + F2{t) , (42) 

where r = -t/iM"^. 

The Regge slope parameter a' can in principle be directly fitted from the knowledge of the 
electromagnetic radii of proton and neutron. In particular, the electric mean squared radii of 
proton and neutron are given by 

4, = .L + ^^. m 

4,„ = rL + ^^, (44) 

where the first term on the rhs is the Dirac radius squared rf, whereas the second term is the Foldy 
term. The Dirac radii are calculated through the integrals of Eqs. (|25|26p for the Rl model, and 
through Eqs. (133|34p for the R2 model. 

In Fig. [H we show the proton and neutron rms radii as the functions of the Regge slope 
a' for both Rl and R2 models. One notes that the neutron rms radius is dominated by the 
Foldy term, which gives r'^ n ~ ' 0-126 fm^. Therefore, a relatively wide range of values a are 
compatible with the neutron data. However for the proton, a rather narrow range of values around 
a = 1.0 — 1.1 GeV~^ are favored. Such value is close to the expectation from the near universal 
Regge slopes for meson trajectories, therefore supporting our Regge type parametrizations. 

In Figs. [21 [21 we show the proton and neutron Sachs electric and magnetic form factors. One 
observes from Figs. [2 [3] that the modified Regge model R2 gives a rather good description of all 
available form factor data for both proton and neutron in the whole t range using the parameter 
for the Regge trajectory a = 1.105 GeV~^, and the following values for the coefficients governing 
the a; — > 1 behavior of the E-type GPDs: ry„ = 1.713 and r]d = 0.566. The 2-parameter version of 
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FIG. 1: Proton and neutron electric mean squared radii p (upper panel) and r'^ „ (lower panel), 
Eqs. (|43|44p . Dotted curves: Regge ansatz according to Eqs. p5|26p : solid curves : modified Regge ansatz 
according to Eqs. p3l34p . Both calculations are shown as function of the Regge slope a . For the quark 
distributions, the MRST02 NNLO parametrization [io'l at scale /i^ = 1 GeV^ was used in the calculations. 
The shaded bands correspond to the experimental values. Note that for the neutron, the Foldy term (term 
proportional to k„ in Eq. (j44|) ) gives f 1; „ = - 0.126 fm^. 



the R2 model gives a description of similar quality if we take a = 1.09 GeV^^ and r/„ = % = 1.34. 

In Figs. [2] and [3l we also show the results of the initial Regge model Rl, with the above value a 
= 1.105 GeV~^ of a'. One sees from Figs. [21 [3] that the Regge model Rl is able to reproduce the 
main trends of both proton and neutron electromagnetic form factor data for —t < 0.5 GeV^. For 
higher values of —t, however, it falls short of the data, since as we discussed, it predicts faster power 
fall-off than that corresponding to the DY relation. The modified Regge model R2 reproduces the 
DY powers for the form factors at large —t, and is able to accurately describe existing data. The 
two additional parameters r/^ and rjd in the R2 model, in particular, allow to describe the decreasing 
ratio of G^/Gj^ with increasing momentum transfer, as follows from the recent JLab polarization 



experiments [38l . l39l . |40| |. Our parametrization leads to a zero for at a momentum transfer of 
— t ~ 8 GeV^, which will be within the range covered by an upcoming JLab experiment 56 1. 

To study the large —t behavior of our GPD parametrizations, it is instructive to plot the Dirac 
and Pauli form factors. In this way, one separates the large —t behavior of both the GPDs H and 
E. In Fig. m we show this large —t behavior for Ff, and for the ratio of F2 /Ff. One observes 
from Fig. |4] that for Ff, the Regge parametrization R2 settles to an approximate ~ 1/t^ power 
behavior around — f ~ 10 GeV^. 

The ratio FfZ-^f ^-Iso discussed within the context of perturbative QCD (pQCD), where the 
asymptotic large-t behavior of the nucleon form factors is dominated by diagrams with two hard 
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FIG. 2: Proton magnetic (upper panel) and electric (middle panel) form factors relative to the dipole form 
Guit) = — as well as the ratio of both form factors (lower panel). The dotted curves correspond 

to the Regge parametrization Rl, with a — 1.105 GeV~^. The solid and dashed curves correspond to two 
fits using the modified Regge parametrization R2. The solid curves are for the 3 parameter fit : a = 1.105 



GeV- 
circles) 



Vu 



1.713 and rjd — 0.566. The dashed curves are for the 2 parameter fit 



a 



1.09 GeV- 



?7d^l.34. Data for the proton magnetic form factor Gj^are from [SlJ (open squares), [32| (open 
(solid stars), ^ (open stars), 35 1 (soHd circles), (sohd squares), according to the recent 
37[. Data for the ratio G^/G^/ are from [38| (solid circles), 39| (open triangles), and 



re-analysis of Ref. 
[iol (solid triangles) 



gluon exchanges [57|, [Sg]. In any model with dimensionless quark- gluon coupling constant, these 
diagrams give Ff ~ 1/t^ [st']. Furthermore, for vector gluons, the quark helicity conservation at the 
gluon vertex and dimensional counting suggest the extra m'^/t suppression for the form factor 



57|,|59|, with m being the quark mass or a nonperturbative parameter coming from the baryon wave 



function corresponding to extra unit of orbital angular momentum [60|. Thus, one should expect 
that F2/F1 ~ 1/t in pQCD. Direct calculation [60], however, shows that the integrals over the 
quark momentum fractions Xi,yj in the pQCD formula contain terms like (p{xi, . . .)(p{yj, . . .)/x?y| 
that diverge even if the nucleon distribution amplitudes (p{xi, . . .), fiuj, ■ ■ ■) linearly vanish at small 
Xi,yj. Strictly speaking, this means that pQCD factorization is not applicable to calculating F^(t) 
even in the asymptotic — t — > 00 limit, the fact well known since the pioneering papers 28|, l59l |. 
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FIG. 3: Neutron magnetic form factor relative to the dipole form (upper panel), and neutron electric form 



The data for the neutron magnetic form factor 
(open triangels), [4j| (solid triangles), [45| (open 



factor (lower panel), with curve conventions as in Fi 

are from [ill (open circles), (solid circles), 
squares), and [46| (solid squares). The data for the neutron electric form factor are from different double 
polarization experiments at MAMI (triangles 47, 48, 49, 5o'|), NIKHEF (sohd square '51]) and JLab (solid 
circles [52j, 53, 54 1). The open squares are the G% extraction from the deuteron quadrupole form factor 
according to the analysis of Ref . [55| . 



The authors of Ref. [gu] substituted the logarithmic divergences by log(— t/Ag^^) factors, and 
obtained pP^^^ /jtP^^^ \og^{—t/K^(jjj)/{—t). This result was found to be in surprisingly good 
agreement with the JLab data. In this connection, we want to emphasize that our results for F2{t) 
and Fi{t) correspond to the Feynman mechanism, i.e., to overlap of soft wave functions. The pQCD 
terms correspond to two iterations of the soft wave functions with hard gluon exchange kernels. 
As is well known, there is ©(a^/vr) suppression for each extra loop of a Feynman diagram in QCD. 
Thus, from our point of view, pQCD terms are ©((a^/vr)^) or, at most, a few per cent corrections 
to the Feynman mechanism contributions to Fi and F2. For this reason, we neglect them in 
our analysis. In our parametrization R2, the good description found for the ratio F2/Ff can be 
directly assigned to the extra suppressing factor of (1 — x)'^ contained in the GPD E{x,t). The 
question, how this suppression is related to the quark orbital angular momentum, deserves further 
investigation. It is interesting to note that the extra (1 — x) factor for £'"'{x) function compared to 
7{^{x) appears in the starting term of the QCD sum rule calculation of these functions 6ll]- Also, 
the dominant x — > 1 perturbative QCD term for the GPD E (given by diagrams) involves two 
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FIG. 4: Proton Dirac form factor (upper panel) multiplied by P and ratio of Pauli to Dirac form factor 
multiplied by —t (lower panel), with curve conventions as in Fig. [21 The data for Ff are from [s^ (solid 
squares). Data for the ratio F^/Fl are from (solid circles), [s^l (open triangles), and 40 1 (solid triangles). 



additional powers in [1 — x) compared with the pQCD expression for the leading a; — > 1 term in 
the GPD H [623. 

Since the GPD E enters the sum rule for the total angular momentum carried by a quark 
of flavor q in the proton as 0] : 

2J'i = J\xx {H^{x, 0, 0) + E'i{x, 0, 0)} , (45) 

our parametrization R2, in which the x — > 1 limit of E is determined from the i^gZ-^f form 
factor ratio, allows to evaluate the above sum rule. The first term in the sum rule of Eq. ()45p 
is already known from the forward parton distributions and is equal to the total fraction of the 
proton momentum carried by a quark of flavor q (q = u, d, s) : 

M| = dxxH^x, 0,0), 

= [ dxx [q^{x) +2q{x)] , (46) 
Jo 

with q{x) the anti-quark distribution. For the 'non-trivial' contribution to the sum rule, arising 
from the second moment of the GPD E, we use our modified Regge parametrization R2 of Eq. (|36p 
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for £'^{x), which, neglecting the antiquark contribution, yields for Eq. 

pi 



2r = + j dxx{l-x)'^^u^{x), (47) 

J^u Jo 

2J'^ = + ^ dxx{l- x^'^ d^{x) , (48) 

2r = M|. (49) 





M| (MRST2002) 


2 J9 (R2 model) 


2 J9 (lattice [63]) 


u 


0.37 


0.58 


0.74 ± 0.12 


d 


0.20 


-0.06 


-0.08 ± 0.08 


s 


0.04 


0.04 




u + d + s 


0.61 


0.56 


0.66 ± 0.14 



TABLE I: Estimate of 2 J' (second column) for the different quark flavors at the scale fi^ — 2 GeV^ 
according to Eqs. (|47ll49p . using the R2 parametrization (with 3 parameters) for the GPD E. For the 
forward parton distributions, the MRST2002 NNLO parametrization 29| is used, yielding the total quark 
momentum contributions (first column). For comparison, the third column shows the quenched lattice 
QCD results of [Sj], extrapolated to the physical pion mass, for 2 J" and 2 J'^. 



In Table HI we show the values of the quark momentum sum rule M| at the scale fi^ = 2 GeV^, 
using the MRST2002 parametrization [29[] for the forward parton distributions. We also show the 
estimate for J", J'^, and J" of Eqs. (fiTlliOl) at the same scale. As was already observed in Ref. [6|], 
based on a Regge model of the type Rl, our estimates lead to a large fraction (63 %) of the 
total angular momentum of the proton carried by the u-quarks and a relatively small contribution 
carried by the d-quarks. As the d-quark intrinsic spin contribution is known to be relatively large 
and negative ( Ad^ ~ —0.25 ), the small total angular momentum contribution J'^ of the ti-quarks 
which follows from our parametrization implies an interesting cancellation between the intrinsic 
spin contribution and the orbital contribution L'^ ( with 2J'^ = Aq+ 2L'? ), which should therefore 
be of size 2L'^ ~ 0.2. For the n-quark on the other hand, the parametrization R2 yields only 
a small value for 2L", as our estimate for 2 J" is quite close to the intrinsic spin contribution 
Auy ~ 0.6. Such a picture is also supported by a recent quenched lattice QCD calculation [63 1 
(see also 6j] for an earlier calculation) for the valence quark contributions to 2 J" and 2J'^. One 
indeed sees from Table H] (third column) that the quenched lattice QCD calculation yields quite 
similar values for 2J" and 2J'^ as our parametrization R2. It remains to be seen however how 
large is the sea quark contribution to the GPD E which can enter the spin sum rule of Eq. (j45p . 
This sea quark contribution is only approximately included (i.e. fermion loop contributions are 
neglected) in the quenched lattice QCD calculations of Ref. 63|]. An exploratory investigation using 
unquenched QCD configurations has been performed in Ref. [65|]. The sea quark contribution is 
also not constrained by the form factor sum rules considered in this paper, which only constrain 
the valence quark distributions. Ongoing measurements of hard exclusive processes, such as deeply 
virtual Compton scattering, provide a means to address this question in the near future. 
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Besides the electromagnetic form factors for proton and neutron, the Regge parametrizations 
discussed in this work can also be used to estimate ^ A transition form factors, provided one 
can relate the N ^ A transition GPDs to the N ^ N ones. First experiments which are sensitive 
to the — > A GPDs have recently been reported 66|]. For the magnetic A'^ ^ A transition form 
factor G^j-(t), it was shown in Ref. [67| that, in the large Nc limit, the relevant — > A GPD can 
be expressed in terms of the isovector GPD E"^ — E'^, yielding the sum rule 



Gl,{t) 



F^-F^it) 



(50) 



where 3.70. Within the large Nc approach used in Ref. [67|], the value G\.j{0) is 

given by G^/j(0) = Ky j [gl, which is about 30% smaller than the experimental number. In our 
calculations, we will therefore use the phenomenological value G^j(O) ~ 3.02 



eg] 




t (GeVO 



FIG. 5: The N ^ A magnetic transition form factor, relative to the dipole form (multiplied by a factor 3). 
Curve conventions as in Fig.[2l The data for G\j are from the compilation of For the JLab data points 
at 2.8 and 4 GeV^, both the analyses of |7o| (upper points) and [iil (lower points) are shown. 



We show our results for using Eq. ([50}) in Fig. [5l It is seen that both the Regge and modified 
Regge parametrizations yield a magnetic A^ ^ A form factor which decreases faster than a dipole, 
in qualitative good agreement with the data. 

The sum rule (|50p was used earlier by P. Stoler 63] , who proposed a model 18|, |l9| in which 
the Gaussian ansatz for GPDs is modified at large —t by terms having a power-law behavior. 



VII. GPDS IN IMPACT PARAMETER SPACE AND POSITIVITY CONSTRAINTS 

The models for GPDs should satisfy many constraints. In fact, such constraints as the re- 
duction of GPDs to usual parton densities in the forward limit and to form factors in the local 
limit, are the key points for the models constructed in this paper. There are more complicated 
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constraints imposed, e.g., by the polynomiality condition which is extremely important for nonzero 
skewness. Since the nonforward parton densities correspond to = 0, they are not affected by 
these constraints. However, they are affected by the positivity conditions which should be taken 
into account both for nonzero and zero skewness parameter. In particular, there exists a relation 
between the ii^-type and H-type GPDs [30|]. Since we are constructing ii^-GPDs from iif-GPDs by 
a simple modification of the x-behavior of -fT by a power of (1 — x), we should check that such a 
modification is consistent with the positivity constraint of Ref. 30l |. 

The most convenient formulation of the positivity constraint relating the E and H GPDs is 
in the impact parameter space. For ^ = 0, the impact parameter versions of GPDs are obtained 
through a Fourier integral in transverse momentum q± : 



(2vr)2 
(27r)2 



-qi) 



(51) 
(52) 



These functions have the physical meaning of measuring the probability to find a quark which 
carries longitudinal momentum fraction x at a transverse position in a nucleon, see Refs. B,[l3]- 
It has been shown 30(] that the GPDs H and E in the impact parameter space satisfy the 
positivity bound : 



2Mn 



\Vk^E\x,h^)\ < Hi{x,hj 



(53) 



Translating the GPD parametrization R2 of Eqs. (j29|35p . into the impact parameter space, we 
obtain : 



E'^ix,h^) 
from which it follows that 

|Vb^ E''{x,h^)\ 



^(1 



j-bxV[-4a {l-x)\nx] 

47r [— a'(l — x) Inx] 



g-bxV[-4a (l-x)lnx] 

Air [— a'(l — x) Inx] 



-bx^/[-4a {l-x)lnx] 



2 



-a 



(1 — x) In X 



2 ' 



(54) 
(55) 

(56) 



Within the R2 parametrization, the positivity bound of Eq. (j53p implies an upper bound on the 
value of I b I I : 



< 



N, 



Mn 



x)'^i 



-a 



[\ — x) Inx 



(57) 



In Fig. [6l we show the GPDs in the impact parameter space for the modified Regge parametriza- 
tion R2 discussed above. The parameters are taken from the best fit to the form factors as discussed 
in the previous section. We see from Fig. [6] that for the n-quark GPDs, the positivity bound of 
Eq. (I53p is satisfied over most of the x-region, considering that the GPDs are vanishingly small for 
values of larger than the nucleon size (corresponding with about 4.4 GeV~^). For the d-quark 
GPDs on the other hand, there is a violation in the present parametrization, which becomes more 
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pronounced at larger values of x and b_|_ , as is shown in Fig. [7] (left panel) . We therefore tried to 
extend the range of validity of the R2 parametrization by finding a fit with a higher value of 77^. 
This can be obtained by imposing the constraint r/„ = r/^. We have shown before that the resulting 
two-parameter fit (a = 1.09 GeV~^, = % = 1.34) gives a nearly as satisfactory description 
of the form factors. It is seen from the right panel in Fig. [7] that this 2-parameter fit extends the 
region in x and b_L for the d-quark where our parametrization satisfies the positivity condition. 

It is clear, that with a somewhat more complicated model, we can easily satisfy the positivity 
constraint. However, given that the violation is rather small, we prefer not to introduce extra 
parameters and to keep the parametrization as simple as possible. 

Furthermore, it is clearly seen from these images that for large values of x, our quark distribu- 
tions are concentrated at small values of b_L, refiecting the distribution of valence quarks in the 
core of the nucleon. On the other hand, at small values of x, the distribution in transverse position 
extends much further out. This expected correlation assures that our model correctly reproduces 
the gross features of the nucleon structure as expressed in terms of the quark distributions. 

VIII. CONCLUSIONS 

Summarizing, we discussed in this work several parametrizations for the i-dependence of the 
nucleon GPDs in view of the recent accurate data for the nucleon electromagnetic form factors in 
the spacelike region. 

Starting from the low —t region, we discussed a Regge model in which the x and t dependence 
of the GPDs are coupled in the form x""'*. This model has only one parameter which physically 
corresponds to the slope a' of the Regge trajectory in the vector EM current channel. This 
parameter is linearly related to the rms radii of F\ and F2 form factors, and it was found that both 
radii are well deccribed by the same universal Regge slope. 

Such a Regge model leads however to faster power fall-off of form factors in the large —t region 
than that expected from the Drell-Yan relation. To conform with this relation and the observed 
power behavior at large —t, we used a modified Regge parametrization that gives slower decrease 
with —t. The modified Regge parametrization displays approximately a l/t^ behavior for Ff{t) 
data in the region -t > 10 GeV^. To describe F^{t), we need to introduce, in addition to a' , two 
parameters that govern the x — > 1 behavior of the GPD E. They were adjusted to give an accurate 
description of the recent polarization data for the ratio Fg/Ff . Since this behavior in our model 
is correlated with the x 1 behavior of the GPD E, it also allows us to evaluate the sum rule 
for the total angular momentum carried by the quarks, which involves the second moment of the 
GPD E. 

For the quark contributions to the nucleon spin, we find an intriguing flavor dependence, in which 
the valence u-quark contributes about two-thirds of the proton's spin (at a low renormalization 
point), which is nearly entirely arising from the u-quarks intrinsic spin contribution. For the valence 
d-quark on the other hand, our parametrization implies a near cancellation between its negative 
intrinsic spin contribution and its orbital angular momentum contribution. Recent quenched lattice 
QCD calculations support this observation. 

It remains to be seen by how much the sea quarks affect this picture. Ongoing measurements of 
hard exclusive processes, such as deeply virtual Compton scattering, are a means to address this 
question. As the GPDs mostly enter in hard exclusive observables through convolution integrals, 
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If, If^, Ff 




012345 012345 
b (GeV"^) b (GeV^) 

FIG. 6: GPDs in impact parameter space for the modified Regge parametrization R2, with parameters 
a = 1.105 GeV~^, rju = 1.713 and rjd = 0.566. Left panels for w-quark; right panels for d-quark. The 
solid (dashed) curves give GPDs H"^ {E'^), respectively. The dotted curves correspond to the function 
|Vb^ -E^l /{2Mn) entering the positivity bound of Eq. (|53l) . 

our parametrization, which builds in the constraint coming from the first moment through the 
nucleon electromagnetic form factors, can be used as a first step to unravel the information on 
GPDs from the observables. The present work also suggests several interesting directions for 
future research. One of them is the extension of this study to quantify the link between the 
nucleon strangeness form factors and the s-quark distributions. Furthermore, the study of the 
chiral corrections (pion mass dependence) to the GPDs will allow to match onto the corresponding 
known chiral behavior of the elastic form factors at small momentum transfer. 
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FIG. 7: Upper bound on bj^ ( Eq. ([57]) ) required by the positivity bound of Eq. ((53l) for two fits using the 
Regge parametrization R2. Left panel is for the 3 parameter fit : a = 1.105 GeV~^, rju — 1.713 and r/d = 
0.566. Right panel is for the 2 parameter fit : a = 1.09 GeV~^, rju = i]d = 1.34. 
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